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Quantifying jet transport properties via large pt hadron production 
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Nuclear modification factor Raa for large pt single hadron is studied in a next-to-leading order 
(NLO) perturbative QCD (pQCD) parton model with medium-modified fragmentation functions 
(mFFs) due to jet quenching in high-energy heavy-ion collisions. The energy loss of the hard 
partons in the QGP is incorporated in the mFFs which utilize two most important parameters to 
characterize the transport properties of the hard parton jets: the jet transport parameter qo and the 
mean free path Ao, both at the initial time tq. A phenomenological study of the experimental data 
for Raa{pt) is performed to constrain the two parameters with simultaneous /A.o.i fits to RHIC 
as well as LHC data. We obtain for energetic quarks go « 1.1±0.2 GeV^/fm and Ao ~ 0.4±0.03 fm 
in central Au + Au collisions at •,/snn = 200 GeV, while go ~ 1.7±0.3 GeV^/fm, and Ao ~ 0.5±0.05 
fm in central Pb + Pb collisions at ^sjvat = 2.76 TeV. Numerical analysis shows that the best fit 
favors a multiple scattering picture for the energetic jets propagating through the bulk medium, 
with a moderate averaged number of gluon emissions. Based on the best constraints for Ao and 
To, the estimated value for the mean-squared transverse momentum broadening is moderate which 
implies that the hard jets go through the medium with small reflection. 

PACS numbers: 12.38.Mh, 24.85.-|-p, 25.75.-q 


I. INTRODUCTION 

A strongly coupled quark gluon plasma (sQGP) con¬ 
sisting of deconfined quarks and gluons may have been 
created in the central region of high-energy nuclear colli¬ 
sions at the BNL Relativistic Heavy Ion Collider (RHIC) 
and the CERN Large Hadron Collider (LHC). One im¬ 
portant evidence for the formation of sQCP from the ex¬ 
periment results are the jet quenching phenomena [Il,i 
that include the strong suppression of single hadron 
spectra , back-to-back dihadron [3, 0 and photon- 
triggered hadron [i,0 as well as jet productions [ll|,Ull 
with large transverse momentum in central A + A colli¬ 
sions as compared to p + p collisions. These observed jet 
quenching patterns in heavy-ion collisions at RHIC/LHC 
can be described well by different theoretical models [13- 
that incorporate parton energy loss induced by mul¬ 
tiple parton scattering and gluon bremsstrahlung as it 
propagates through the dense matter. 

A widely used parameter controlling the parton en¬ 
ergy loss is the jet transport parameter q |25l| . or the 
mean-squared transverse momentum broadening per unit 
length for a jet in a strong interacting medium, which 
is also related to the gluon distribution density of the 
medium [^. and therefore characterizes the medium 
property as probed by an energetic jet. To estimate the 
jet transport parameter q intense theoretical studies have 
been made, such as with the weakly-coupled QCD [13- 
113], the strong-coupled AdS/CFT [lo,[3l|, and the lattice 
calculations (321 l33l | . Recently, a phenomenological inves¬ 
tigation had been carried out to extract the initial value 
of jet transport coefficient qo at initial time tq, which 
gives go ~ 1-2 ± 0.3 CeV^/fm in Au-|-Au collisions at 
y/sNN = 200 CeV and go ~ 1.9±0.7 CeV^/fm in Pb-fPb 
collisions at = 2.76 TeV for a given quark with ini¬ 

tial energy of 10 CeV [33| . 


In this paper we will extract the initial jet transport 
parameter go and the initial mean free path Ao at ini¬ 
tial time Tq on the bulk medium evolution by compar¬ 
ing the experimental data at RHIC/LHC with numerical 
simulations of of single hadron yields with large px in a 
next-to-leading order (NLO) pQCD parton model, where 
the EPS09 parametrization set of NLO nuclear parton 
distribution functions (nPDFs) has been used to take 
into account of possible initial-state cold nuclear mat¬ 
ter effects, and a phenomenological model [13, O for 
medium-modified fragmentation functions calculated in 
leading-order (LO) at twist-4 in the high-twist approach 
of jet quenching [35l - l37| | has been utilized to incorporate 
parton energy loss. The evolution of bulk medium used 
in the study for parton propagation was given by a 3 -I- 
1 dimensional ideal hydrodynamic model [38l l39j which 
is constrained by experimental data on hadron spectra. 
From calculations with the two independent inputs for 
the parameters and simultaneous y^/d.o.f fits to the 
RHIC and the LHC data, we obtain that for energetic 
quarks go ~ l.I ± 0.2 CeV^/fm and Aq ~ 0.4 ± 0.03 
fm in central Au -I- Au collisions at y/SNN = 200 CeV, 
while go ~ 1.7 ± 0.3 CeV^/fm, and Aq ~ 0.5 ± 0.05 fm 
in central Pb + Pb collisions at ^/snn = 2.76 TeV. This 
simultaneous and separate constraint of the two initial 
values should give a precise and quantitative description 
for jet quenching to probe the medium properties. For 
a parton jet propagating through the bulk medium, the 
average transverse momentum broadening squared (g^) 
depends on the transport parameter as well as the mean 
free path, (g|,) = gA. Our numerical results show that 
the mean transverse momentum broadening squared of 
energetic partons for one scattering at initial time tq in 
the center of the fireball at LHC is about 2 times of that 
at RHIC. 

The rest of the paper is organized as follows. We first 
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give a brief overview of the NLO pQCD parton model for 
single inclusive hadron spectra and a phenomenological 
model for medium-modified fragmentation functions in 
Sec. El Then the numerical calculations for phenomeno¬ 
logical studies of the experimental data on single hadron 
suppression and extraction of the jet transport parame¬ 
ter and the mean free path are carried out in Sec ■ ED 
We present some discussions in Sec. m and finally sum¬ 
marize our study in Sec. El 


II. NLO PQCD PARTON MODEL AND 
MODIFIED FRAGMENTATION FUNCTIONS 


We will utilize the pQCD parton model at NLO for 
the initial jet production spectra which has been applied 
to large px hadron production in high energy hadron- 
hadron reactions with great successes [i^. In the model 
the differential cross section of hadron yields has been ex¬ 
pressed as a convolution of NLO parton-parton scatter¬ 
ing cross sections, parton distribution functions (PDFs) 
in nucleons and parton fragmentation functions (FFs), 
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where da{ab —>■ cd)/dt denotes the leading-order (LO) 
elementary parton scattering cross sections at a^. The 
NLO contributions in O(a^) involve both 2 —>■ 3 tree 
level processes and one loop virtual corrections to 2 —>■ 2 
tree processes. Processes at 2 —3 tree level inlcude 
qq -t qqg, qq -t qqg, qq -t ggg, qg -t qgg, qg -t qqq, 

gg —> qqg, gg —>■ ggg, etc, which include soft and collinear 
contributions. A standard MS renormalization scheme 
is applied to control ultraviolet divergence in one loop 
virtual corrections to 2 —2 tree processes. More de¬ 
tailed discussions on calculations at NLO could be found 
in . In this paper the numerical calculations are car¬ 
ried out with a NLO Monte Carlo program where 
two cut-off parameters, 5s and 5c, are employed to iso¬ 
late the collinear and soft divergences in the squared ma¬ 
trix elements of the 2 —>■ 3 processes. The regions with 
the divergences are integrated over in n-dimension phase 
space and the results are added with the squared matrix 
elements of the 2 —>■ 2 processes. This gives a set of two- 
body and three-body weights depending on 5s and 5c. 
But the dependence will be eliminated after the weights 
are combined in the calculation of physical observables, 
and the hnal numerical results are insensitive to the cut¬ 
off parameters [40j| . 

We employ the same factorized form for the inclusive 
large px particle production cross section in nucleus- 
nucleus collisions, which can be computed as a convo¬ 
lution of nuclear thickness functions, the nuclear par- 
ton distribution functions (nPDFs), elementary parton- 
parton scattering cross sections and effective medium- 
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at fixed impact parameter b in the transverse plane of 
the beam direction. In Eq. the average over the az¬ 
imuthal angle of the initial fast parton is implicitly im¬ 
plied. The nuclear thickness function t(r) is calculated 
with the Woods-Saxon distribution function for nucle¬ 
ons in a nucleus and has been normalized by requiring 
f d'^rtAi'x) = A. The nuclear parton distributions per 
nucleon (nPDFs) fa/A{xa,y'^,'x) can be parameterized 
as the production of the parton distributions inside free 
nucleons fa/N{x, fj?) and the nuclear shadowing factor 

Sa/A{x,y?',Y), 
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where Z denotes the charge and A is the mass num¬ 
ber of the nucleus. In the numerical simulations we 
use the CTEQ6M parametrization for nucleon par- 
ton distributions f„/v(a;, m^), and EPS09 parametriza¬ 
tion of nPDEs [^. Since the parton-parton scattering 
cross sections are computed up to NLO, the CTEQ6M 
parametrization and EPS09 parametrization are both 
used at NLO. Eor simplicity, we only use the central- 
fit set of EPS09 parametrization in following numerical 
calculations. 

An energetic parton jet produced in the hard scatter¬ 
ing may suffer multiple scattering with thermal partons 
in the QGP created in nucleus-nucleus collisions. The 
jet-medium scattering and medium induced gluon radi¬ 
ation should give rise to new contributions to parton 
fragmentation functions (FFs) in vacuum and thus leads 
to medium-modified fragmention functions, which may 
evolve with the scale Q [13,113, [13, m, Ii3| , in a similar 
way like the DGLAP evolution in vacuum. If the par- 
ton jet travels a distance L inside the medium with the 
inelastic scattering mean free path A, the probability for 
the jet scattering n times to the medium can be assumed 
to obey Poisson distribution [13, Therefore, effect 
of parton energy loss in the dense QGD medium can be 
calculated in the high-twist approach of jet quenching 
and the effective medium-modified parton fragmentation 
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functions (mFFs) can be given by BE Em, 




^ nO /,, 2 \ 
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where z' = pr/ipTc — ^E^) is the rescaled momen¬ 
tum fraction of the hadron from the fragmentation of 
the quenched parton which has the initial transverse 
momentum pTc and loses energy AFlc during its prop¬ 
agation inside the hot medium, z'g = {L/X)pt/ A.Ec is 
the rescaled momentum fraction of the hadron from the 
fragmentation of a radiated gluon with initial energy 
A_Ec/(j), and Zc = PtIptc is the momentum fraction 
for jet fragmentation in vacuum, (j) times scattering 
will provide (j) gluon emissions, so there is a factor 
(j) in the fragmentation contribution of emitted gluon 
in above equation. As shown in Ref. Bl, the above 
mFFs satisfy the momentum sum rule by construction, 
Eh f zDh/c(zc,fJ.^,AEc) = 1. 

The weight factor exp(—(j)) is the probability for 
those partons escaping the medium without suffering any 
inelastic scattering, and the weight factor 1 — exp(—(j)) 
is the probability for partons encountering at least one 
inelastic scattering. The rescaled fraction in Eq. o is 
got by energy shifting due to energy loss AEc. For a 
given jet, the energy loss AEc and the scattering number 
(j) both depend on the local medium density in the jet 
trajectory and characterize the medium properties. This 
approximative approach for medium-modified fragmen¬ 
tation function reproduces the main effect of medium- 
induced radiation [l4| , and is therefore similar to another 
approximative approach in references 1^ where the 
modified fragmentation function is concluded as a con¬ 
volution of the vacuum fragmentation function and the 
probability for a given jet to be quenched to a final jet 
inside the medium. We also note the difference between 
our approach and the ones used in references [3113 ■ In 
higher-twist formalism , one considers twist-4 pro¬ 

cesses of the splitting of a highly virtual parton (/i >> 
Aqcd) in QCD medium and evaluates the contribution of 
medium-induced gluon radiation, which gives rise to the 
effectively modified parton fragmentation functions and 
their corresponding (medium-modified) QCD evolution 
equations with respect to the hard scale p. This is dif¬ 
ferent from that in references (dsl . 113 where the medium 
contribution is computed at the medium scale. Also in 
references [13113 n convolution with a Poisson probabil¬ 
ity for multiple emissions is used for the energy shift to 
take into account the fluctuation of energy loss. 

In a high-twist approach the total energy loss AE = 
AEc{E,h,r) is related to the jet transport parameter 
via, 



xqF{y)4:Sm^{xLP~''y /2), (5) 


as shown in Ref. [53, in which one can refer for more 
details. y~ denotes a jet place in its trajectory, and is 
the same as normal time r used in following calcula¬ 
tions. In the above expression, the LPM effect for the 
induced gluon emission originates from the destructive 
interference of two kinds of processes, i.e. the soft-hard 
and hard-hard processes, which become identical to each 
other and lead to a cancellation of their contributions 
when the transverse momentum of the radiated gluon is 
small and the formation time of radiated gluon is rather 
large [ssl-ls^- This is exactly the same as the LPM effect 
in the double scattering in the GLV opacity expansion 
formalism [^. In the case of multiple soft scattering, 
as discussed in references [13113 lE^-Bl the LPM effect 
is caused by similar interferences and plays a dominant 
role when the coherent emission of a single gluon in a 
multiple scattering process is considered and leads to a 
suppression of the energy loss as compared to the addi¬ 
tive contribution of n{= (j)) independent scattering of 
one gluon radiation. 

We emphasize that in our model, the parton matrix 
elements are calculated at NLO, and EPS nuclear PDFs 
at NLO are utilized to include the initial-state cold nu¬ 
clear matter effects, whereas the parton energy loss due 
to the final-state hot medium effect is given by a LO 
derivation of parton energy loss resulting from effective 
medium-modified fragmentation functions in twist-4 at 
the high-twist expansion approach. Recently several the¬ 
oretical attempts [23 - 1^ . Issl . [13 have been made to cal¬ 
culate momentum broadening and parton energy loss due 
to multiple scattering in medium beyond leading order, 
which may be incorporated phenomenologically to a com¬ 
plete NLO calculations of particle or jet productions in 
high-energy nuclear collisions. It is also noticed that re¬ 
cently theoretical investigation of color decoherence has 
been developed for jets resol ving and energy redistribut¬ 
ing in the QCD medium 

The jet transport parameter for a gluon is 9/4 times 
of a quark, which is assumed to be proportional to the 
local parton densiW in a dynamical evolving medium and 
expressed as [13 [^ [Mj , 


g(r, r) = qo 


Pg(r,r + nr) pt^Uf, 

Pg (7-0,0) po 


( 6 ) 


for a parton produced at a transverse position r at an 
initial time tq and traveling along the direction n. go 
denotes the jet transport parameter at the center of the 
bulk medium at the initial time tq . Pg is the gluon density 
at a given temperature r(T, r), and in numerical calcu¬ 
lations we assume pg oc (1 — f)T^ for the medium as an 
ideal gas. As introduced in Ref. [H, jH, [^ , the fraction 
/(r, r) of the hadronic phase at any given time and local 
position is given by 


f 0 

/(T,r) = .1 0 - 1 


if T>170 MeV, 
if T=170 MeV, 
if T<170 MeV 


( 7 ) 
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FIG. 1: (Color online) (a) The hadron cross sections at mid- 
rapidity in p+p and central Au + Au collisions at ^/snn = 200 
GeV. (b)(c) The corresponding nuclear modification factor 
with different values of the jet transport parameter qo and 
the mean free path Aq. The data are from |63l - [65l| . 


FIG. 2: (Color online) (a) The hadron cross sections at mid¬ 
rapidity in p+p and central Pb+Pb collisions at ^sjvjv = 2.76 
TeV. (b)(c) The corresponding nuclear modification factor 
with different values of the jet transport parameter qo and 
the mean free path Ao- The data are from [1,0|. 


where we consider the mixed phase contribution and ne¬ 
glect the pure hadron phase contribution. During the 
mixed phase at T = 170 MeV the hadron phase frac¬ 
tion will be / = 0 ^ 1 while the QGP phase faction 
will be 1 — /. With the time evolution during the mixed 
phase the fraction value / will increase from 0 to 1. In 
the following numerical calculations the time-dependent 
fraction f is given by simulations of the hydrodynamic 
model |38l.l39l|. is the four momentum of the jet and 
is the four flow velocity in the collision frame. The aver¬ 
age number of scatterings along the parton propagating 
path is given by, 


(L/A) 


J_ /■°“ Pg(T,r-FnT) 
AoX„ ^ Pg(ro,0) ’ 
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where Aq is the mean free path at the initial time Tq, 
and for a quark jet it is 9/4 times of that for a gluon 
jet. The parameter Aq as well as go will be indepen¬ 
dently inputted in the following numerical calculations 
and simultaneously constrained by experiment data. 

The fragmentation function in vacuum ZI°y^(zc, in 
Eq. dl]) and (|3|) is given by the AKK parametrization [fi^ . 


Here we use the NLO AKK FFs parametrization. For in¬ 
clusive hadron production given by Eq. o and there 
are three independent scales: the renormalization scale 
fj-rem the factorization scale pfact and the fragmentation 

scale Pfrag- choOSe pren — f^fact — l^frag — ^•‘^PT la 

our numerical analysis for both p+p and A+A collisions. 


III. EXTRACTING PARAMETER VALUES IN 

SIMULTANEOUS xVd-o.f FITS TO Raa DATA 

In the model for the effective medium-modified frag¬ 
mentation functions (mFFs) disucssed in Sec. IIIIl infor¬ 
mation on the space-time evolution of the local temper¬ 
ature and flow velocity in the bulk medium along the jet 
propagation path should be provided. In our simulation 
we will utilize a (3-1-1) dimensional ideal hydrodynamic 
model m, 11^ to obtain the space-time evolution of the 
bulk matter created in central nucleus-nucleus collisions. 

With a given space-time profile of the gluon density, 
one can then utilize the preceding effective mFFs to ob¬ 
tain the high px hadron spectra. In actual calculations 
for the spectra or cross section at fixed values of the 
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hadron transverse momentum px in Eq. 0 and m, 
the factorization and renormalization scales are all cho¬ 
sen as pf = fiR = 1.2 pt- Shown in Fig. [T] (a) are the 
hadron cross sections \n p + p and 0-5% Au Au col¬ 
lisions at y/SNN = 200 GeV with given parameter val¬ 
ues go = 1-1 GeVVfm and Aq = 0.4 fm, as compared 
to PHENIX data [63l - [6^ . The theoretical cross section 
ioi A + A collisions is scaled as {dN aa/ dydPpT)/TAA{b) 
where TaaQ}) = / d'^rtA{r)tA{r — b). We choose 6=2 
fm for 0-5% Au + Au and b = 2.1 fm for 0-5% Pb + Pb 
collisions. Shown in Fig. [5] (a) are for 0-5% Pb + Pb colli¬ 
sions with given parameter values go = 1-7 GeV^/fm and 
Ao = 0.5 fm, as compared to CMS and ALICE data [1,0. 
It is observed that the theoretical results with chosen pa¬ 
rameters of go and Ao could describe the experimental 
data at the RHIC and the LHC very well. 

Shown in (b)(c) panels of Fig. [T] and Fig. [5] are the 
suppression factor or nuclear modification factor. 


Raa 


dNAA/dyd'^PT 
TAA{b)dapp/dyd'^PT' 


(9) 


To compare theoretical results with the RHIC and LHC 
data, we may fix one parameter of go or Aq, and then 
choose different values for another parameter. For given 
Ao the nuclear modification factor decreases with the in¬ 
creasing of go, while for given go the nuclear modification 
factor increases with the increasing of Ao. 

For the two parameters go and Aq we fix one parameter 
and constrain another by x^/d.o.f fitting to data for the 
nuclear suppression factor. The y^/d.o.f is defined as 
follows, 


xVd-O.f 




(Vth - Ve^pf 






( 10 ) 


where Vth stands for the theoretical value, V^xp denotes 
the experimental value, gives the quadratic sum 

over all types of errors that one chosen point has, and N 
the number of data points selected. 

In numerical calculations the jet transport parameter 
go and the mean free path Aq are two independent inputs. 
We choose for a quark jet go = 0.1 - 3.0 GeV^/fm and Ao 
= 0.1 - 1.0 fm, while for a gluon jet the values are cor¬ 
respondingly multiplied by 9/4 for go and 4/9 for Ao be¬ 
cause of different color factors in gluon-gluon and quark- 
gluon interacting vertex. From simultaneous x^/d.o.f fits 
to experimental data at RHIC and LHC shown in Fig.|31 
one can extract values of the jet transport parameter go 
and the mean free path Aq at the center of the most cen¬ 
tral A-\- A collisions with the given initial time tq = 0.6 
fm. For a energetic quark jet, the best fits to the com¬ 
bined PHENIX data [63l - l^ give go = 1.1±0.30 GeV^/fm 
and Ao = 0.4±0.03 fm in 0-5% central Au-\- Au collisions 
at y/SNN = 200 GeV, while the best fits to the com¬ 
bined ALICE 0 and CMS 0 data lead to go ~ 1.7 ± 0.3 
GeV^/fm, and Ao ~ 0-5 ± 0.05 in 0-5% central Pb-|-Pb 
collisions at ^snn = 2.76 TeV. 


vVd.o.f. for R„„ in 0-5% Au+Au at 200 GeV 

K aA 




xVd.o.f. for R„ in 0-5% Pb+Pb at 2.75 TeV 



(GeV^/fm) 



FIG. 3: (Color online) The x^/d.o.f. as a function of the 
initial quark jet transport parameter go and the initial mean 
free path Aq. Upp er panel is from fitting to the combined 
PHENIX data [634^ on Raa{pt) for 7 r ° with pr ~ 5 - 20 
GeV at mid-rapidity in 0 — 5% central Au + Au collisions at 
^ysNN = 200 GeV. Lower panel is from fitting to the combined 
ALICE 0 and CMS i data on Raa{pt) for charged hadrons 
with Pt = 10 - 100 GeV at mid-rapidity in 0 — 5% central 
Pb + Pb collisions at ^/snn = 2.76 TeV. 


IV. DISCUSSIONS 


In general, the jet transport parameter should depend 
on the scale as shown by recent studies on the renormal¬ 
ization of the jet transport parameter [bbl - lh^ . In our 
numerical simulations for the initial values for go and Ao 
given by Eq. m and Eq. (|S]), we assume they are con¬ 
stants for different jet transverse momentums as a rea¬ 
sonable approximation for phenomenological studies at 
the RHIC and the LHC. Note that the best fits shown in 
Fig. [3] are obtained for hadrons with px = 5-20 GeV at 
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the RHIC while pt = 10-100 GeV at the LHC, so what 
we constrain for and Aq may be understood as the 
averaged values for jets with different transverse momen¬ 
tum. Roughly speaking, what we constrain for go and Aq 
should depend on much greater transverse momentum at 
the LHC than at the RHIC because of the much wider 
kinematical region of jet pT at the LHC. 

In the formulation for the medium-modified fragmen¬ 
tation functions in Eq. o, the final-state medium effect 
of jet quenching is controlled both by the total energy loss 
AE oc go in Eq. ([5]) and the multiple scattering number 
(j) oc ^ in Eq. Therefore the suppression factor 
Raa is quantified by the two independent parameters go 
and Ao. The trend for the simultaneous x^/d.o.f fits in 
Fig. [3] shows that an increasing go must associate with 
an increasing Ao to fit well the data at both the RHIC 
and the LHC. In fact, a larger Aq gives a smaller scat¬ 
tering number, and then a larger go is needed to release 
greater energy loss per scattering in order to describe 
experimental data. This implicit relation between the 
two parameters is consistent with theoretical estimates 
for the jet transport parameter and the mean free path 
which are related each other via the local temperature in 
a weakly-coupled QCD medium [2^ . 

Of interest are the two different limits as demonstrated 
in Eig. [21 the single scattering limit with large Ao, e.g. 
Aq = 0.7 fm at RHIC, and the infinite number scatter¬ 
ing limit with very small Ao, e.g Ao = O.I fm at LHC. 
The numerical simulations for the simultaneous y^/d.o.f 
fits at both RHIC and LHC show that in the single scat¬ 
tering limit the suppression factor Raa is insensitive to 
go and sensitive to Ao, whereas in the infinite number 
scattering limit Raa is sensitive to go and insensitive to 
Aq. According to an assumption [il, for parton scat¬ 
tering obeying a Poisson distribution, the probability for 
those partons escaping the system without suffering any 
inelastic scattering is exp(—(j)), while the probability 
for partons encountering at least one inelastic scatter¬ 
ing gives I — exp(—(j)). One can see these two weight 
factors in the medium-modified fragmentation function 
in Eq. (|3]). In the infinite number scattering limit with 
small Ao, exp(—(j)) is very small with large j, the first 
term of Eq. (|4|) with dependence on go will dominate the 
total fragmentation contribution, so Raa is sensitive to 
go and less insensitive to Ao. On the other hand, in sin¬ 
gle scattering limit with large Ao, the second term of Eq. 

gives the dominant contribution, thus Raa is insensi¬ 
tive to go and more sensitive to Ao. We observe that our 
best fits for go and Ao are found in the region between 
the single scattering limit and the infinite number scat¬ 
tering limit due to competing effect between the energy 
loss per scattering quantified by go and the scattering 
number quantified by Ao, which implies that the data fa¬ 
vor a regime of mean-free-paths that suggests multiple 
scattering in the medium. 

Erom Fig. 121 we can extract go range of values for ener¬ 
getic quarks as constrained by the measured suppression 
factors of single hadron spectra at RHIC and LHC as 


200GeV AuAu^/ 0-5% y=0.4 fm 
2.76TeV PbPb^h= 0-5% y=0.5 fm 
Jet Collaboration 


0.35 


0.40 

T (GeV) 


FIG. 4: (Color online) The scaled jet transport parameter 
q/T^ for an initial quark jet at the center of the most cen¬ 
tral A+A collisions at an initial time ro = 0.6 fm extracted 
by comparing the theoretical simulations with experimental 
data at both RHIC and LHC. The dashed boxes indicate the 
corresponding results of Jet Collaboration [j^ ]. 


follows: 
go Ri 


I.I ± 0.2 ir 

I 7± 0 3 


T=373 MeV, 
T=473 MeV, 


at the highest temperatures reached in the most central 
Au-|-Au collisions at RHIC and Pb-|-Pb collisions at LHC. 
As shown in Fig.|2]for the scaled jet transport parameter 
q/T\ our result falls within the range of go/T q ener¬ 

getic quarks extracted from experimental data on Raa 
by Jet Collaboration though it is considerably smaller 
than that given by a strong-coupled AdS/CFT calcula¬ 
tion [s^, El as well as a lattice calculation [s^, [13 ■ 

In addition, from Fig.|2]one can extract Ao range of val¬ 
ues for energetic quarks as constrained by the measured 
suppression factors of single hadron spectra at RHIC and 
LHC as: 


, / 0.4 ±0.03 „ T=373 MeV, 

^ \ 0.5 ± 0.05 T=473 MeV. 

In a theoretical estimate [13 for the mean free path, 
1/Ag = pa = 2>as{Q^)T for which the elastic cross sec¬ 
tion tr is used at leading order and the density p is for 
an ideal gas. One can introduce K factor to account for 
higher order correction and the more realistic interaction 
among the medium particles 0[nl, 

^ = SKa,iQ^)T. (II) 


Considering A, = |Ag, and assuming the scale = ET 
for a hard parton with energy E traversing a hot QCD 
medium with temperature T, we find with K = 2.5 - 
4.0 Xq given by above equation is equal to our best fit 
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^ 0.6 - 
C3 


<q;> = q>. 


200GeV AuAu -> n 0-5% = 0.6 fm 

2.76TeV PbPb ^ h= 0-5% t. = 0.6 fm 


T (GeV) 


FIG. 5: (Color online) The averaged number of gluon emis¬ 
sions from a propagating quark as a function of the vertex 
place of the hard scattering. The created quarks propagate 
along -l-x direction to escape off the fire ball. 


for the mean free path at the highest temperatures in 
both RHIC and LHC. The K factor is bigger than what 
would be naturally expected, which might be caused by 
LO a as well as the ideal gas density p in the theoretical 
evaluation of 1/Ag = pa. Higher order correction for the 
cross section may provide a factor of ~ 2, while other 
effects such as corrections due to the difference between 
the real dynamics of the QGP and the simple picture of 
an ideal gas may account for the remaining enhancement 
to K. For instance, a strongly interacting QGP may give 
a larger cross section than a weakly-coupled QGP. Thus 
the comparison of the model simulation with the data 
seems imply that the hot QGD medium at the RHIG and 
LHG is more likely a strongly interacting medium, which 
is surely mo del-dependent and further validations from 
other observables will be needed for a robust conclusion. 
It is noted that in numerical estimates we use parton 
energy = 8 - 25 GeV for pr = 5 - 20 GeV hadron 
production at RHIG while i? = 15 - 120 GeV for pT = 
10- 100 GeV hadron production at LHG, and therefore 
the running coupling as(Q^} is appreciably smaller at 
LHG than at RHIG. 

The phenomenological model given by Eq. iD assumes 
that one scattering will induce one gluon emission from 
the propagating parton, so for a given propagating par- 
ton the total scattering number equals to the total num¬ 
ber of gluon emissions from this parton. Recent theo¬ 
retical calculation gives the averaged number of gluon 
emissions (Ng) from a propagating parton in HT-BW 
approach within the high-twist framework of parton en¬ 
ergy loss [lllzl- In the HT-BW model, the medium- 
modified FFs are given by numerically solving a set of 
modified DGLAP evolution equations within the high- 
twist approach with an initial condition given by Poisson 
convolution of multiple gluon radiations, which has been 


FIG. 6: (Golor online) The temperature dependence of the 
average transverse momentum broadening squared for ener¬ 
getic quarks for one scattering at the initial time tq in the 
center of the fireball. 


shown to give the best agreement with data for the 
nuclear modification factor Raa hr high-energy heavy-ion 
collisions. Especially, the averaged number of gluon emis¬ 
sions (Ng) from a propagating parton is given in the 
study for modified DGLAP evolution equations, which 
can be compared with our extracted number of medium- 
induced emissions (j). Shown in Fig. [5]is the compari¬ 
son for gluon emission number between our model (solid 
and dash curves) and the HT-BW approach (dot and 
dotted-dash curves denoted as “HT-BW”), where the ini¬ 
tial quark jets are produced in the point (x, y = 0) of x 
axis and propagate along -\-x direction in the transverse 
plane to escape off the fire ball. The quark transverse 
momentum is for example chosen as = 10 GeV for 
RHIG and 30 GeV for LHG in central A-l-A collsions. 
Our results for the averaged number of gluon emissions 
are consistent with the HT-BW method, and justify the 
validity of the model as shown in Ref. [13, [l3 ■ 

According to definitions for the jet transport parame¬ 
ter and the mean free path [^ . 

Q = pJdqT-^qT: ( 12 ) 



one can estimate the average transverse momentum 
broadening squared, 

(9t) = (14) 

Then our best fitting values for qo and Aq give, 

/ 0.44±0.11 2 T=373 MeV, 

[qj,) _ qA Ri I Q gg ^ T=473 MeV, 






for energetic quarks with one scattering at the initial time 
tq in the center of the fireball, as shown in Fig.[51 Our nu¬ 
merical results show that for energetic parton jets scatter¬ 
ing inside the bulk medium at the highest temperature, 
the average transverse momentum broadening squared at 
LHC is about twice of that at RHIC. Compared to initial 
parton jet energy, the broadening is moderate, which im¬ 
plies that the jet may traverse through the medium with 
small reflection and justifies the eikonal approximation 
usually used in jet quenching calculations. 

As we stated before the AKK FFs in vacuum are 
used in our numerical simulations. It is noted that a 
recent theoretical study [t^ has confronted seven sets 
of NLO FF parameterizations with inclusive charged- 
particle spectra in p-l-p collisions at the LHC and iden¬ 
tified that most of the theoretical predictions including 
AKK08 tend to overpredict the measured cross sections 
by up to a factor of two due to the too-hard gluon-to- 
hadron FFs. In this paper we focus on the medium 
properties demonstrated by the nuclear modification fac¬ 
tor Raa which is a ratio of spectra between A + A and 
p + p collisions and therefore is expected not to be very 
sensitive to the choice of FFs parametrization as well as 
the scale. We have redone our simulations with Kretzer 
parametrization of FFs [t^, which show the extracted go 
(Ao) is less (larger) about 10-20% by using AKK08 FFs 
than by using Kretzer FFs. 


V. SUMMARY 

We have used the NLO pQCD parton model with ef¬ 
fective modified fragmentation functions due to radia¬ 


tive parton energy loss to study single hadron spectra in 
high-energy heavy-ion collisions at both RHIC and LHC. 
The energy loss of the hard partons is incorporated in 
the modified fragmentation functions which utilize two 
most important parameters to characterize the proper¬ 
ties of the bulk medium, the jet transport parameter go 
and the mean free path Aq both at the initial time tq. 
We perform the phenomenological study of the experi¬ 
mental data for Raa{pt) to constrain the two parame¬ 
ters with simultaneous x^/d.o.f fits to RHIC as well as 
LHC data, and obtain for energetic quarks go ~ 1.1 ± 0.2 
GeV^/fm and Ao ~ 0.4 ± 0.03 fm in central Au -I- Au 
collisions at ^snn = 200 GeV, while go ~ 1.7 ± 0.3 
GeV^/fm, and Ao ~ 0.5 ± 0.05 fm in central Pb + Pb col¬ 
lisions at ^ysNN = 2.76 TeV. Numerical analysis shows 
that the best fit falls between the single scattering limit 
and multiple scattering limit for the energetic jets prop¬ 
agating through the bulk medium. These results indi¬ 
cate that the average transverse momentum broadening 
squared (g^) = gA of energetic partons for one scatter¬ 
ing at initial time tq in the center of the fireball at LHG, 
{Qt)LHC ~ 0.85 GeV^, which is twice of {q^)RHic ~ 0.44 
GeV2 at RHIC. 
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